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ABSTRACT 



We explore the prospects for searches at the LHC for sparticle decays that violate r lepton 
number, in the light of neutrino oscillation data and the seesaw model for neutrino masses and 
mixing. We analyse the theoretical and phenomenological conditions required for tau flavour 
violation to be observable in \2 — > X + decays, for cosmologically interesting values of 

the relic neutralino LSP density. We study the relevant supersymmetric parameter space in 
the context of the Constrained Minimal Supersymmetric Extension of the Standard Model 
(CMSSM) and in SU(5) extensions of the theory. We pay particular attention to the possible 
signals from hadronic tau decays, that we analyse using PYTHIA event simulation. We find 
that a signal for r flavour-violating X2 decays may be observable if the branching ratio exceeds 
about 10%. This may be compatible with the existing upper limit on r — * fij decays if there 
is mixing between right-handed sleptons, as could be induced in non-minimal SU(5) GUTs. 
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1 Introduction 



Data from both atmospheric [1] and solar [2] neutrinos have confirmed the existence of neu- 
trino oscillations with near-maximal — v T mixing (Super-Kamiokande) and large v e — > 
mixing (SNO). These observations would also imply violation of the corresponding charged- 
lepton numbers, which would be enhanced in supersymmetric theories and might be observable 
in low-energy experiments. In fact, charged- lepton- number violating processes could occur at 
embarrassingly large rates if the soft supersymmetry-breaking masses of the squarks and slep- 
tons were not universal. For this reason, it is often assumed that these masses are equal at 
the grand-unification scale, as in the constrained minimal supersymmetric extension of the 
Standard Model (CMSSM). 

Even within the minimal supersymmetric version of the seesaw model for neutrino masses, 
renormalization of the soft supersymmetry-breaking slepton masses would occur while running 
from the GUT scale to the heavy neutrino mass scales. This would be induced by the Dirac 
Yukawa couplings of the neutrinos [3], since these cannot, in general, be diagonalised simul- 
taneously with the charged-lepton and the slepton mass matrices. This scenario provides a 
minimal amount of charged-lepton-flavour violation, which could be further enhanced by GUT 
interactions and/or non-universal slepton masses at the GUT scale. 

Within this framework, many signatures for charged-lepton-flavour violation have been con- 
sidered [4, 5], including ji — > decays, ji — e conversions, r — > /ry and r — * decays. In 
view of the (near-) maximal mixing observed amongst the corresponding neutrino species, one 
expects these decays to be relatively large when the soft supersymmetry-breaking masses Mi/ 2 
and/or m are relatively small. Other charged-lepton-flavour violating possibilities that have 
been considered are the decays X2 — ¥ X + e ± /i =F [6, 7], and X2 —> X + Ai ± T =F [8, 9], where x 
is the lightest neutralino, assumed here to be the lightest supersymmetric particle (LSP), and 
X2 is the second-lightest neutralino. It has been argued that these decays might have a rate 
observable at the LHC for certain choices of the CMSSM parameters [8, 9]. These decays would 
provide search possibilities that are complementary to searches for the flavour-violating decays 
of charged leptons, since they may be relevant particularly for regions of the supersymmetric 
parameter space where rare charged lepton decays and conversions are suppressed. 

The answer to the question which decay mode offers better detection prospects at the LHC 
or at a future linear collider depends on the details of both the theoretical model and the 
experiment. The decay X2 X + Ai ± T =F has an experimental signature that is less distinctive 
than X2 X + e ± Ai =F - However, it may have certain theoretical advantages over the latter mode. 
This is because v T — mixing is known to be essentially maximal, and the feedthrough into 
the charged-lepton sector is potentially enhanced by larger Dirac Yukawa couplings and/or 
lighter singlet-neutrino masses, if neutrino masses exhibit the expected hierarchical pattern. 
These comments imply that the two lepton-flavour-violating modes are complementary, and 
both have to be studied. In order to assess the observability of X2 X + r± / i=F ) the detailed 
simulation of signal and background events seems unavoidable. 

In previous work [9] , we observed that the branching ratio for X2 —> X + T ± Ai =F decay is enhanced 
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when m X2 > m fl > m x , where fi is the lighter stau slepton. This occurs in a wedge of 
the (M1/2, mo) parameter plane in the CMSSM that is complementary to that explored by 
t — > /ry. The region of CMSSM parameter space where this enhancement occurs includes the 
region where x~ t coannihilation suppresses the relic density Q x , keeping it within the range 
Q x h 2 ~ 0.1 preferred by astrophysics and cosmology, even if Mi/ 2 is comparatively large. The 
interest of this coannihilation region is supported by experimental constraints on the CMSSM, 
such as m h and b — > 57 decay, which disfavour low values of My 2 . 

In the current paper, we revisit the X 2 — * X + t ± /U t decay mode and extend previous analyses 
in the following directions: 

• Theoretical framework. In [9] we considered generic seesaw mixing within the framework of 
the CMMSM. Here, we will study how the ratio T(x2 X + r± + ^ T ) l^(X2 — > X + r± + r=F ) is 
modified under different theoretical assumptions, and what information a possible signal may 
provide in this respect. Specifically, in our previous analysis we considered only the dominant 
effects due to large mixing in the 2-3 sector of the charged-slepton mass matrix. Within this 
framework, the slepton mixing arises essentially only from the LL sector, and large LFV effects 
may be observed mostly at low tan/?, where the smuon and stau are almost degenerate. As 
tan/3 increases, the LFV width decreases if flavor is only violated in the LL sector, since the 
lightest stau becomes mostly right-handed. However, in a GUT-inspired model the running 
of couplings from the Planck to the GUT scale may introduce significant corrections to the 
right-handed slepton masses, giving in principle rise to the possibility of enhanced rates at 
large tan f3 as well. We address these considerations in a more elaborate study of the complete 
mixing effects. 

• Detailed study of the supersymmetric spectrum and parameter space: We pay particular atten- 
tion to regions that lead to large values of T( X 2 X + r± + r=F ) through on-shell slepton produc- 
tion, (BR( X 2 - XtV) = Yti [BR(X2 - l^)BR(h - r X ) + BR( X 2 - hr)BR(h - l x X )]), 
while satisfying all phenomenological and cosmological (relic density) constraints. The charac- 
teristic parameter region for the signal in the r channel to be optimal is defined by the following: 

(i) m X2 > m f > m x ; 

(ii) we also assume that one of the mass differences in (i) is > m T and the other > m M , mf > m x , 
so that the //, r and f are all on-shell; 

(iii) Moderate values of m x (phase space and luminosity considerations). 

These conditions are obeyed in significant fractions of the stau coannihilation region. 

• PYTHIA event simulation: We use PYTHIA to simulate the hadronic decays of rs produced in 
the dilepton decay of the second- lightest neutralino, X2 ~* M ~* [<*>]• in the study of the 
flavor- violating dilepton signal (r ± /i =F ), the second lepton is tagged as a muon with a probability 
equal to the branching ratio assumed for flavor-violating decays. 

After first setting the theoretical scene, our procedure in this paper is a bottom-up one, namely: 

(i) Knowing what branching ratio would be required for an observable signal at the LHC, 
through an event simulation performed using PYTHIA; 

(ii) We study the theoretical frameworks that may satisfy the conditions for observability, 
together with the supersymmetric parameter space requirements. In doing so, we review look 
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at models with mixing arising mainly in the left-handed slepton sector, which could not produce 
an observable signal, and then discuss the possibility of mixing in the right-handed sector, which 
becomes possible in non-minimal GUTs, where effects due to renormalization above M GUT 
schemes may become relevant. 

The structure of our paper is as follows. Different theoretical frameworks are reviewed in Section 
2, the supersymmetric parameter space is explored in Section 3, our simulation is described in 
Section 4 and its results in Section 5. Our conclusions are summarized in Section 6. 



2 Possible Sources of Charged-Lepton Flavour Violation 

2.1 Renormalization below Mqxjt 

In the minimal supersymmetric extension of the seesaw mechanism for generating neutrino 
masses with three heavy singlet-neutrino states N iy the leptonic sector of the superpotential is: 

W = N^(Y u )ijLjH 2 - E^Y^jLjH, + ^N^MijN^ + //# 2 #i , (1) 

where Y v is the neutrino Dirac Yukawa coupling matrix, M.^ is the Majorana mass matrix 
for the Ni, the Lj and Hj are lepton and Higgs doublets, and the are singlet charged- 
lepton supermultiplets. The superpotential of the effective low-energy theory, obtained after 
the decoupling of heavy neutrinos is [10] 

W eff = LiH 2 (yJ (M ) 1 yJ\ LjH 2 — El ( Y e ) ^ LjE\. (2) 

\ / ij 

In the basis where the charged leptons and the heavy neutrino mass matrices are diagonal, one 
finds 

M V = Y? (M D y l Y v v 2 sin 2 (3, (3) 
where v = 174 GeV and tan/3 = v 2 jv\. 

In the context of the CMSSM, the soft supersymmetry-breaking masses of the charged and 
neutral sleptons are assumed to be universal at the GUT scale, with a common value m . In the 
leading-logarithmic approximation, the non-universal renormalization of the soft supersymmetry- 
breaking scalar masses is given by 

implying that the corresponding corrections to left-handed slepton masses are given by Vd, the 
Dirac neutrino mixing matrix in the basis where the ci-quark and charged-lepton masses are 
diagonal. In this approach, non-universality in the soft supersymmetry-breaking left-slepton 
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masses is much larger than that in the right-slepton masses, particularly when the trilinear 
soft supersymmetry-breaking parameter A = 0. Consequently, within the CMSSM, to a 
good approximation, the renormalization of the soft supersymmetry-breaking parameters at low 
energies can be understood in terms of the dominant non-universality in the third-generation 
left-slepton mass: mg M = diag(mo, m 2 ,, x x wig), where a typical value of the non-universality 
factor is x ~ 0.9. Based on the above, in [9] we assumed that there is an off-diagonal tl — Jil 
mixing term in the soft mass-squared matrix 



Am, 



Oll 



[1 -x)m 



2 sin(20) 



(5) 



where is the mixing angle between the second and third generation in the charged-lepton 
Yukawa matrix. This mixing leads to lepton-flavour violation ~ sin 2 (20). 

In our current work, we go beyond the above approximation by including the complete mixing 
effects in interesting regions of the parameter space. For the structures of the mixing matrices, 
as well as the heavy and light neutrino hierarchies, one must appeal to a specific GUT model. 
Here we consider textures obtained by combining SU(5) with a U(l) family symmetry. Requir- 
ing that fields in the same GUT multiplets have the same flavour charge leads straightforwardly 
to fermion mass matrices of the following forms [11, 12] 



M u oc 



e 6 e 5 e 3 
e 5 e 4 e 2 



M down oc 



e 4 e 3 e 3 
e 3 e 2 e 2 



Me OC 



e 4 e 3 e 



e 3 e 2 



e 3 e 2 



(6) 



with e ~ 0.2 as the preferred value. The mass structure and the mixings of the neutrinos 
are more complicated, because of the heavy Majorana masses of the right-handed components 
arising from terms of the form z/rz/rE, where S is an SU(3)x SU(2) xU(l)-invariant Higgs field 
with Iw = and a non-zero flavour U(l) charge. A working example was given in [13]: 



M, 



rr oc 
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(7) 



Appropriate choices of the unspecified U(l) charges, such as n\ = 2,n 2 = — 1,^3 — l,a — —1 
(where the n» are the U(l) charges of the right-handed neutrinos, and a is the U(l) charge of 
the field S) lead to interesting phenomenology. 

We use an indicative choice of coefficients for these SU(5) textures given in [13], which is 
summarised in Table 1 below. 
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Parameters in an SU(5) model with large tan/3 


Charged leptons 


a\ 2 = 0.6, a e 13 = 0.9, a e 22 = 1.2, a e 23 = -0.5e i7r / 3 

°31 = 0-7; = 0.6, 033 = 0.4 


rrf D 


^ = 1.3,^ = -1.3,a^2 = 0.7 
= l.Se-/ 5 , ^2 = 0.7,^3 = 0.5 


M RR 


0^2 = 1, a 33 = 1-8 



Table 1: Choice of coefficients that reproduce the fermion data for an SU(5) model with large tan/3. 
Coefficients not listed in the table are set to unity. 



2.2 Renormalization above Mqut 

One must also keep in mind that the GUT scale may lie significantly below the scale Mc ra v at 
which gravitational effects can no longer be neglected 1 . In this case, the renormalization of 
couplings at scales between Mcrav and Mqut may induce significant flavour-violating effects, 
particularly in the right-handed slepton mixing, which is suppressed in minimal schemes such 
as that described in the previous subsection. The simplest such example is provided by the 
minimal supersymmetric SU(5) GUT, where the superpotential contains terms of the form 
EUH (with H being a colour-triplet Higgs field that is expected to have a mass ~ Mqut)- 
This gives rise to one-loop diagrams that renormalize the right-handed slepton masses. In the 
leading-logarithmic approximation, these corrections take the form [14]: 

(mjh - -^ A «3^^'(3^ + ^)log^, (8) 

for % ^ j, where Vu denotes the mixing matrix in the corresponding couplings in the basis 
where the w-quark and charged-lepton masses are diagonal. Similarly, the complete leading- 
logarithmic renormalization of the A e terms is given by 

- -^ a o(A ei U^^A^log^ + A ei U-V^A^log^ 

+3A e ,^^A^ 3 log^). (9) 

For the structure of the mixing matrices Vu,d, one has to go to a specific GUT model. Within 
the minimal SU(5) GUT, the ci-quark mass matrix is the transpose of the charged-lepton mass 
matrix, Vd is simply the unit matrix, and Vu is related to the familiar CKM matrix. This 
mechanism therefore provides relatively little mixing. 



1 Mcrav might be identified with either the Planck mass Mp = 1.2 x 10 19 GeV or some lower string unification 
scale M st ring ~ 10 18 GeV. 
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2.3 Non-Minimal GUT Effects 



In the minimal SU(5) case, the renormalization would be too small to generate observable 
mixing for universal initial conditions. However, minimal SU(5) also predicts the unsuccessful 
relations m s = and m d = m e . These can be modified by non-renormalizable terms in the 
effective superpotential, such as the fourth-order term 10 — 24 — 5 — H, which make different 
contributions to the rf-quark and charged-lepton mass matrices [15], such as: 

A(10 - 5 - H) + X'(H - 10 - 24 - 5) -> Xv(dd c + e c e) + X'vV(2dd c - 3e c e). 

In this case, in the basis where rrid is diagonal, one has m e = — hX'vV , where the matrix 
of couplings A' is non-diagonal, in general. The diagonalization of = V e Rm e V^ L then gives 
= V e Ft{m d — bX'vV)V^ L . Similarly [16], the colour-triplet-induced e c u c mixing may receive 
large corrections for the first two generations. Parametrizing the non-renormalizable correction 
to this mixing by V u r, the RGE-induced right-slepton mixing is not given by Vqkm as in the 
minimal model, but by the product Vr = VckmVuRi where V u r is a potential extra source of 
right-handed charged-lepton mixing, that can be significantly larger than what is expected in 
the minimal scheme [15] 2 . 

3 Discussion of the Supersymmetric Parameter Space 
3.1 First Considerations 

In order to identify representative points in the supersymmetric parameter space for which 
detailed simulations are to be run, we first recall a point considered previously in the litera- 
ture [8] 3 . In this work, the following CMSSM point has been studied: 

tanf3 = 10, m = lOOGeV, M 1/2 = 300GeV, A = 300GeV. (10) 

This selection of parameters is displayed as point A in Table 2. 

We consider first the effects of lavor-mixing entries on the slepton mass matrices that are 
introduced to mimic the non-diagonal terms induced in the M\ L sector by a generic 'seesaw' 
machanism: 

K) 2 3= 5 'K) 22 - < n > 

2 These non-renormalizable corrections also change the forms of the fcrmion mass matrices, and hence the 
predictions of this type of flavour-texture model within minimal SU(5). Thus, these corrections would also affect 
the renormalization between the GUT and heavy-neutrino mass scales. Such effects would have supplementary 
effects on the the left-handed slepton mixing, but the detailed study of those effects goes beyond the scope of 
this paper. 

3 This point is now excluded by the LEP bound on the Higgs mass. Nominally, nih > 114.4 GeV, but we 
consider supersymmetric points with as low as 111 GeV to be acceptable, so as to make a suitable allowance 
for uncertainties in the theoretical calculation of m^. 
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Point 


Modeltype 


m 


M 1/2 


tan (3 


A 


^events 


"int 


Lint 


A 


CMSSM 


100 


300 


10 


300 


757K 


25.3 pb 


30 fb- 1 


B 


SU(5) 


40 


450 


35 


40 


730 K 


2.44 pb 


300 fb" 1 


C 


CMSSM 


220 


500 


35 


220 


536 K 


1.79 pb 


300 fb- 1 



Table 2: Parameters of the two reference points A and B (masses in GeV). We also quote the 
numbers of events simulated, the LHC cross sections and the assumed sample luminosities. Point A 
is a CMSSM model with universal soft supersymmetry-breaking terms at the GUT scale. Point B is 
a model with universality assumed at a scale 2 • 10 17 GeV; for comparison with this point, we also 
present point C, a set of CMSSM parameters that leads to a similar sparticle spectrum and satisfies 
all the cosmological and phenomenological bounds. In all cases, we work with jj, > 0. 



Point 


M- g 


Mu L 




X.?. 


M fl 


M~o 

Xi 


Mr 


Mr 


M h 


A 


720 


664 


669 


216 


150 


118 


155 


232 


110 


B 


1095 


1025 


1024 


366 


207 


194 


286 


371 


117 


C 


1154 


1074 


1078 


388 


219 


206 


290 


405 


116 



Table 3: Relevant sparticle masses (in GeV) for the reference points defined in Table 2. 



including all the contributing on-shell sfermion exchange diagrams, as given in [17]: 

3 

BR( X 2 - XtV) = E [BR(X2 - kv)BR(li - t X ) + BR( X 2 - hr)BR(h - » x j\ • (12) 



i=l 



We present in the left panel of Fig. 1 the ratio of the flavor- violating decay width T(x2 — > 
X + + fi T ) to the CMSSM flavor-conserving decay width T(x2 — * X + r± + rT ) f° r 8 = 0- We 
find similar values to those in [8] , where only the dominant f i exchange diagram was included. 
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As shown later, in order to have significant LFV signals, we need T(x2 X J i~ T +/ i=F )/r(X: 
X + + t t ) ~ 0.1. We see in the left panel of Fig. 1 that this ratio corresponds to 5 ~ 0.25. 
However, a potential problem for such points arises from BR(r — > /ry), which tends to become 
too large. This is shown in the right panel of Fig. 1, in which typical values are much larger 
than the estimated value of BR(r — > 10 -9 in [8]. We recall that the current experimental 

upper bound is < 6.8 • lO -8 [18]. 

When a similar choice of CMSSM parameters is considered in the framework of an SU(5) 
model with seesaw neutrinos like the one described in Sections 3 and 4 we find that by setting 
Mr = 3 ■ 10 14 GeV we obtain a value of 5 ~ 0.05 leading to a prediction BR(r — > /ry) ~ 3 • lO -8 . 
This would be acceptable for r — > /ry, but would not lead to observable X2 ~~ ^ X + T± + ^ T 
decay. Due to the strong bound imposed by r — > ^7, it is very difficult to obtain reasonable 
values of 5 using only the LL mixing found in seesaw models. 

However, mixing in the RR sector can enhance the decays so that they might be detectable 
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5 



Figure 1: In the left panel, flavor-conserving and -violating dilepton branching ratios are calculated 
for point A, for comparison with Fig. 1 of [8]. The corresponding expectations for r — > ^7 decay are 
shown in the right panel as a function of 5. 

at the LHC without a large increase in BR(r — > ^7), since RR mixing enters only in the 
subdominant one-loop neutralino-exchange diagram for this process. It is considerably smaller 
than the dominant chargino-exchange diagram for this process, which is sensitive to LL mixing. 
However, as we have already discussed, significant mixing in the RR sector of the slepton mass 
matrix cannot be obtained in minimal SU(5) with the conventional seesaw mechanism, but may 
be obtained in more general GUT scenarios 4 

In order to see how a branching ratio for % 2 — > X + r± + decay that could be observable at 
the LHC might be induced by non-diagonal entries in the Mrr sector of the charged-slepton 
mass matrix, without violating the r — > /I'y bound, we parametrize as follows the LFV entries 
in the charged-slepton mass matrices: 



We illustrate in Fig. 2 the different dependences of the branching ratios for X2 — * X + T + M T 
and r — > ^7 decays on the magnitudes of the off-diagonal LL and RR LFV mixing entries. 

The upper plots in Fig. 2 display the dependences of the xi ~ ¥ X + r± + £* T decays of interest, 
and the lower plots display the dependences of the branching ratio for r — > /ry decay. The left 
plots show the dependences on 5ll for certain discrete choices of 5rr, and the roles of LL and 
RR mixing are reversed in the right plots. We see again that if RR mixing were negligible a 10% 
ratio of the LFV X2 —> X + r± + ^ T decays relative to the flavour- conserving ones would require 
large LL mixing with 5ll ~ 0.25, which would inevitably imply a violation of the BR(r — > ^7) 
bound. On the other hand, we see that a value of 5rr ~ 0.02 — 0.03 would lead to a similar 
ratio for % 2 — > x + r ± + yU T decay, whatever the value of 5ll, whilst the LFV radiative tau decay 

4 Large mixing in the RR sector was also found to be necessary for flavour-violating effects to be observable 
via non-universality in leptonic kaon decays [15]. 




(13) 
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could remain within the acceptable experimental range if 5ll were below about 0.03. We also 
observe that the size of 5rr required to obtain the desired ratios increases as tan/3 increases, 
requiring physics beyond the minimal SU(5) model, as discussed in Section 2.3. 





A 



Ci -7 

S3 in r- 
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-// Srr= ° 


E\p. bound 






1 0.1 


.2 0.3 
8 LL 


0.4 




Figure 2: Branching ratios for point A for \2 —> X + r± + decay (upper plots) and r — > fij (lower 
plots) as functions of 5ll (left) and 8rr (right) for certain discrete choices of Srr (left) and 5ll 
(right), for the same CMSSM inputs as in Fig. 1. 



3.2 Phenomenological and Cosmological Constraints 

The parameter point A discussed in the previous section predicts a neutralino relic density 
that exceeds the WMAP upper bound, and also predicts m h < 111 GeV, below even the 
more conservative LEP bound on the Higss mass allowing for theoretical uncertainties. In this 
subsection we specify a model that satisfies these phenomenological constraints, and may also 
incorporate the non-minimal SU(5) GUT model with a seesaw mechanism described earlier. 

We assume universal soft supersymmetry breaking at a scale Mx = 2 ■ 10 17 GeV> Mgut, in 
which case we find a minimum value of tan f3 ~ 31 for fi > 0, above which f — x coannihilations 
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Figure 3: Cosmologically-favored areas (green) in the (Mi/ 2 ,mo) plane for tan/3 = 35 and Aq = tuq, 
assuming SU(5) unification. In the left panel we assume universality at Mx = Mqut, whereas in 
the right panel we assume universality at Mx = 2 • 10 17 GeV. The red areas are excluded because 
m x > mf. We also display the contours for = 111,114 GeV (black solid and thin solid) and 
BR{b — > 57) • 10 4 < 2.15,2.85 (blue dashed and thin dashed). 

are sufficient to bring Q x h 2 below the upper WMAP bound [19]. Furthermore, we find that a 
minimum value tan/3 ~ 35 is needed to ensure m/, > 114 GeV, though we recall that we assign 
a theoretical uncertainty of ~ 3 GeV to the theoretical calculation of m^. The GUT Yukawa 
coupling relation Y T = Y b cannot be achieved for /1 > [20] due to the large supersymmetric 
threshold corrections to m&. However, lepton mixing effects may modify the GUT relation: 

Y b (GUT) = Y T {GUT){\ - x), (14) 

where x is a parameter that accounts for a sizeable 2-3 generation mixing in the charged-lepton 
Yukawa coupling matrix [21, 22, 23]. 

We display in Fig. 3 two (Mi/2,mo) planes for the same choice tan/3 = 35. In the left panel, 
we assume universality at Mx = Mqut, whereas in the right panel we show the changes in 
the allowed parameter space when Mx = 2 ■ 10 17 GeV. In the latter case, we find a region 
of parameter space in which the WMAP bound on the cold dark matter density is respected 
simultaneously with the conservative bound > 114 GeV. In both cases, we fix rrib(Mz) = 
2.92 GeV, which corresponds to the evolution of the central value of the MS value m^mb) = 
4.25 GeV with a s (M z ) = 0.172. With this choice of m b , we find a value x = 0.37 in (14). 

In order to take into account the cosmological and phenomenological considerations discussed 
above, we now study point B in Table 2, whose sparticle spectrum is tabulated in Table 3. 
Point A is taken from [8] and is included so as to facilitate comparisons. As we shall see, point 
B allows for an observable number of events with the luminosity expected at the LHC. 

We display in Fig. 4 for point B the same sets of branching ratios as were shown in Fig. 2 
for point A. We see in the top panels that in this case we also need large non-diagonal entries 
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Figure 4: Branching ratios, as in Fig. 2, for the point B described in the text. 

in the slepton mass matrix in order to achieve a branching ratio for x 2 Xi' r± / i=F that is of 
interest for the LHC, e.g., 5 RR ~ 0.15 for 5 LL = or 5 LL ~ 0.35 for 5 RR = 0. We also see in 
the lower panels that r — > ^7 is again very restrictive on the size of 5 l l, imposing a maximum 
value ~ 0.03. We see in the bottom-right panel that S RR ~ 0.15 is allowed for Sll = 5 . The 
choice Mx = 2 • 10 17 GeV opens up the possibility of generating such a value of S RR though 
RGE-induced mixing above M GUT . If we assume a generic pattern of non- universal M RR GUT 
entries of the form: 

I m 2 \ 

M 2 RR (GUT)= m 2 t-m 2 , (15) 
\ e-m„ m| / 

a value e ~ 0.25 is needed to generate the value 5 RR ~ 0.15 at the electroweak scale. 

5 If one chose a model with (nominally) the lower value rrih = 111 GeV, one could allow lighter masses than 
for point B, however the limits from BR(r — > ^7) on 5rr and Sll would be more severe. 



12 



4 Signals and Backgrounds at the Event Generator Level 



This section is devoted to considerations about the detection of LFV in neutralino decays at 
the LHC, at the event generator level. The signal of tau flavour violation that we consider is an 
excess of dilepton pairs of the type fir h over er h pairs, where signifies a hadronic jet produced 
by hadronic tau decay. The hadronic decay channel for the tau lepton provides a good SUSY 
signature because: (i) tau leptons decay 65 % of the time to hadrons, and (ii) hadronic jets 
from tau decays are narrower and have lower particle multiplicity than conventional QCD jets, 
and they can be reconstructed from the tracking and calorimeter information. Although their 
charges cannot be determined reliably, tau leptons decaying hadronically can be detected more 
clearly than taus decaying to lighter leptons because leptons arising from tau decays have two 
neutrinos, and so are more difficult to reconstruct than leptonic tau decays and separate from 
background processes. 

We present, analyse and compare results from two supersymmetric models with different spar- 
ticle spectra, one formulated within the CMSSM and another incorporating additional GUT- 
inspired physics, as summarized in Tables 2 and 3. The spectra for these points were calculated 
using ISA JET 7 . 78 [24] and then interfaced into PYTHIA 6.418 [25] , which was used to generate 
the cross sections, the ensuing QCD parton showers and the hadronic tau decays. 



4.1 General Event Cuts 

There are some general and distinctive features of R-conserving sparticle production at the 
LHC. First, the dominant production mechanism is expected to be the pair production of 
massive squarks and gluinos, which subsequently decay into hadrons, resulting in several hard 
jets. On the other hand, the neutralino LSP escapes detection and thus we expect large missing 
transverse energy (MET) in SUSY events. These features suggest using the following definition 
for the effective mass of the process: M e ff =$ t +Pt,i + Pt,2 + Pt,3 +Pt,4, where p tti stand for the 
transverse momenta of the four hardest jets. The distribution in M e ff would exhibit a peak 
around the scales of the squark and gluino masses. 

Keeping these features in mind, we apply the following cuts on the events generated by PYTHIA, 
for both the SUSY signal and the Standard Model background. 



(i) N jets > 4 with pr tl > 100 GeV and p r ,2,3,4 > 50 GeV, 

(ii) t > 0.2M eff and p t >fi t min , 

(Hi) M eff =^ t + pr tl + p Tt2 + p Tt3 + p T>3 + p TA > M™f GeV; 

The values of fi t mm and M™f were fixed independently for each set of SUSY parameters, in 
order to optimize the suppression of the Standard Model background in each case. 
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4.2 Calorimeter and Object Definition 

Only calorimetric smearing and segmentation have been incorporated, since a full detector 
simulation is out of the scope of the present work. A uniform segmentation A0 = Ar] = 0.1 is 
assumed and the energy smearing is done according to the following detector parametrization: 

ECAL ~ 10% /VE + 1%, | 77 |< 3, 
HCAL ~ 50%/VE + 3%,\r] \<3, 
FCAL ~ 100%/ VE + 7%, | 

We use the following object definitions: 

• Jets are identified with the help of the PYCELL subroutine of PYTHIA, with R = 0.4 used 
for the jet cone size, due to the high jet multiplicity in SUSY cascade events. We take 
cells with more than 1 GeV as possible jet triggers and require a sum of Et > lOGeV in 
the calorimeter for the jet to be finally accepted. 

• The total missing transverse energy, E T , is defined as the vector sum of the deposits in 
the calorimeter cells defined above. 

• Leptons are required to be central in rapidity, with | rj |< 2.5, and we use the following 
isolation demand: no more than 10 GeV of transverse energy should be present in a cone 
size of R = 0.2 around the lepton direction. 

• Hadronic tau decays were selected by requiring > 20 GeV and | rj |< 2.5, and a 
'matching' jet with p T T > 0.8 p 3 ^ 1 and | rj |< 2.5. The jet was tagged as a tau in the two 
following cases: (i) if the jet-tau distance is no more than R = 0.4 in the rj — plane, 
or (ii) if the jet-tau distance R > 0.4, the jet could also be accepted as a tau, with a 
probability 

P = 1 - (0.971j4 /2 - 49) §(1 " £t) (16) 
where e T = 0.7 is the selected efficiency. 

• Since ISA JET docs not allow for LFV decays, flavor- violating dileptons (r ± /x T ) were sim- 
ulated by counting events with two taus, with at least one of them decaying hadronically; 
the second tau was then tagged as a muon with a probability equal to the assumed LFV 
branching ratio (10 %). 

4.3 Standard Model Background 

Standard Model processes that could in principle contribute to the background for the signa- 
ture we study here are ti, WW, QCD jets, Z jets and W jets. The dominant backgrounds 
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are expected to be ti and WW, and only these are considered here. In practice, no event 
coming from the other processes would satisfy the above-mentioned kinematical requirements 
when simulated using PYTHIA, as used here, and a more extensive study of Standard Model 
backgrounds is out of the scope of this paper. Such a study would require the use of specialized 
codes for the evaluation of NLO cross sections and of the exact matrix elements for the parton 
showers and hadronic decays. The background samples and cross sections used in this work are 
summarized in Table 4 6 . 



Process 


10 fb- 1 


100 fb" 1 


a 


ti 


250K 


1M 


486 pb" 1 


WW 


250K 


1M 


70 pb- 1 



Table 4: Standard Model background samples for reference luminosities 10 fb 1 and 100 fb 1 , together 
with their LHC cross sections calculated using PYTHIA 



5 Results for LFV at the LHC 

In this Section, we analyse the sets of parameters described in Table 2, presenting the results 
of our event simulations. We recall that point A has been studied previously in [8], and 
it is included here solely as a consistency check. Point B is based on the SU(5) RGEs in the 
framework discussed in previous Sections, demanding compatibility with phenomenological and 
cosmological bounds. We also recall that point C of Table 2 yields a low-energy spectrum and 
results that are very similar to those of point B. 

Plots corresponding to point A (usually placed at the tops of the figure arrays), show the 
numbers of events normalized to a reference luminosity of 10 fb -1 . Plots corresponding to 
point B (usually placed at the bottoms of the figure arrays) are for an integrated luminosity of 
100 fb" 1 . 

5.1 Choices of Cuts 

In Fig. 5, we show the M e ff and Missing Transverse Energy (MET) distributions of events 
simulated for point A (top) and point B (bottom). 

Plots (a) show the distributions calculated including all the events for both the Standard Model 
and SUSY, i.e., without applying any of the cuts in (16). We see that the SUSY events exceed 

6 Because the simulations are quite time-consuming, rather than generate the total number of events corre- 
sponding to the reference luminosity, we have rescaled the total number of events by an overall factor which, in 
the case of 100 fb -1 , is roughly 2.5 times that for 10 fb -1 . 
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Figure 5: M e ff distributions for the crude events (plots a) and after applying cuts (i) and (ii) (plots 
b), and the Missing Transverse Energy (MET) distributions after applying the same cuts (plots c). 
The plots on the top correspond to simulations for point A, where E™ m = 0; those on the bottom to 
simulations for point B with E^ nm = 320 GeV. The signal distributions are shown as (red) solid lines, 
and the Standard Model backgrounds by the (blue) shaded histograms. 



the Standard Model ones at large values of M e ff > 1200 GeV for point A, whilst the distribution 
for point B shows larger numbers of Standard Model events than SUSY events over the entire 
M e ff range. 

In plots (b), after imposing cuts (i) and (ii) (with = for point A and $™ in = 320 GeV 
for point B, the SUSY events exceed the SM events at « 800 GeV and ~ 1400 GeV for points 
A and C respectively. We use these values for M™J in cut (iii) in the different cases. 

Plots (c) show the MET distribution after applying the above cuts. The SUSY event rates 
exceed those of the Standard Model at energies above ~ 240 GeV for point A and w 480 GeV 
for point B. Thus, no cut in MET is needed for point A (it would reduce the Standard Model 
background but also the SUSY signal, as can be seen below for the invariant mass distributions, 
which we denote by IMD). However, a cut in MET is required for point B, reducing the total 
number of accepted SUSY events for point B relative to point A. 
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5.2 Invariant Mass Distributions 



In Fig. 6 we show (left) dilepton IMDs for light leptons and (right) hadronically-decaying r 
leptons after imposing the cuts explained in (16), as fixed in the previous subsection 5.1. In 
the case of point A (top plots), both the opposite-sign same-flavour (OSSF) light dileptons 
(left) and r ± — r T distributions (right), both shown as red solid lines, exhibit breaks at the 
characteristic end-points 7 , which occur at the invariant mass 



M max 



(M? — Mj)(Mf — M? ) 

, — ^ l — = 98.3 GeV. (17) 

N i 



On the other hand, the opposite-sign different-flavor (OSDF) (left) and the same-sign t^t^ 
distributions (right) do not display such a kinematical constraint. For point B (lower plots) 
the OSSF end-point does not exist in the light lepton case, and the peak at M z is completelly 
obscured by the Standard Model background. In the OS r — r case, the end-point is at 
108.3 GeV. 



5.3 Lepton-Flavour- Violating Events 

The left plots of Fig. 7 show the visible mass distributions for pairs of hadronically-decaying 
taus. Excesses of OS pairs (red solid lines) over the SS pairs (green dashed lines) ending 
approximately at the kinematic end-points can be seen. The structures of the end-points are 
not as clean as in Fig. 6, because of the energy carried away by neutrinos. In the right plots, 
an excess of OS I t t^ pairs over the SS pairs can be seen. The distributions of OS LFV /it 
pairs (that cannot be distinguished experimentally from I t t^) are also shown for the sake of 
comparison. In view of the difficulty in distinguishing experimentally the LFV /i — r signal 
pairs from the Standard Model background, we simulate and plot the excess of fi — r pairs over 
e — r pairs. These should be identical in the Standard Model on average, so any excess of // — r 
pairs beyond statistical fluctuations would be a signal of LFV. 

In fig. 8 (left) the excesses of /i T T^ over e T r^ pairs for points A and B are shown. The 
observable numbers, N l J^, of fJ, T r^ LFV pairs are obtained by summing the counts in the 
subtracted /U T r/f — e T r^ distributions in the interval of M; r masses between 30 and 110 GeV. 
We obtain 

Point A : Nj£ = 470 ± 39 (12 a) 

Point B : N l j£ = 308 ±30 (10 a) (18) 

where we quote only the statistical errors for the signal samples. If we estimate an efficiency 
of 70 % for the jet-tau matching, a lower number of (fJ>Th)if v pairs would be obtained: 

7 Howcver, the solid curve in the top right panel of Fig. 6 is not observable, since the energy carried off by 
the v T would escape detection. 
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Figure 6: (Left) Dilepton invariant mass distributions for light leptons: opposite-sign same-flavor 
(OSSF) (red solid line), opposite-sign different-flavor (OSDF) (green dashed line) and the OSSF Stan- 
dard Model background (shaded). (Right) Hadronically- decaying ditau mass distributions (including 
the v T four-momenta) with OS (red solid) and SS (green dashed), and the Standard Model background 
for OS ditau pairs (shaded). The top plots are for case A, and the bottom plots are for case B. 



Point A : N l / T v h = 355 ± 34 (10 a) 

Point B : N 1 ^ = 236 ± 27 (9 a) (19) 

In fig. 9, the excesses of LFV fiT^ pairs over ei"^ are shown in comparison with the subtraction 
of er h pairs from \iT h pairs obtained from the simulation of the Standard Model background, for 
points A (left) and B (right). The non-zero subtracted signals for the Standard Model result 
from the statistical fluctuations in our simulation. The comparison shows that the statistical 
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Figure 7: Left: visible ThTh mass distributions with OS (red solid lines), SS (green dashed lines) and 
the OS Standard Model backgrounds (shaded). Right: visible Wh mass distributions with OS (red solid 
lines), SS (green dashed lines), and OS Standard Model backgrounds (shaded). The LFV OS [xr^ pairs 
are also shown (pink- dot- dashed). As usual, the top plots are for case A and the bottom plots are for 
case B. 

significance of the LFV signal is quite high in both cases, as it can be separated well from the 
Standard Model statistical fluctuation: we obtain S/B ~ 2.6 for point B, whilst the ratio is 
better for point A, namely S/B ~ 4.65. 



5.4 Results for Varying tuq at Fixed Mi/ 2 

So far, we have restricted our attention essentially to CMSSM points lying in the coannihilation 
strip where the f\ is only slightly heavier than the neutralino LSP, which restricts the kinematics 
of the r hadronic jet. Points with the same value of Mxji but larger values of m may also be 
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Figure 8: Left: comparisons between the visible mass distributions for r^fj, (red solid lines), r^e (green 
dashed lines) and Standard Model ThfJ* pairs (shaded). The LFV fiTh pairs have been added to the fATh 
distribution. Right: comparison between the visible mass distributions for LFV [it^ (red solid lines), 
the sign- subtracted distribution l^rj^ — l^T h (green dashed lines), and the sign- subtracted Standard 
Model backgrounds (shaded). 

allowed under certain circumstances, e.g., in direct-channel resonance regions that may appear 
at large tan f3 in the CMSSM or at lower tan (3 in models with non-universal Higgs masses, or 
if the gravitino is the LSP. Points with larger mo have different kinematics, so here we study 
the consequences for the observability of LFV in such models. We fix M1/2 = 500 GeV and 
tan/? = 35, and vary uiq and (slightly) A . We choose four reference points, whose parameters 
are listed in Table 5, and the resulting sparticle mass spectra in Table 6. 

We display in Fig. 10 the numbers of LFV h^t^ pairs for each of these points, calculated in the 
same way as in previous subsection, for the points described in Tables 5 and 6. We also display 
the statistical error bars. We see that the LFV signal continues to be observable up to the 
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Figure 9: The signal for excess fj,Th LFV pairs (red solid lines) and subtracted \xrh — er^ Standard 
Model backgrounds (shaded) for points A (left) and B (right). As no Standard Model contributions 
should survive after subtraction, on average, the shaded signals are due to statistical fluctuations. 



Point 


m 


Mi /2 


tan/3 


A 


N events 


0~int 


Lint 


1 


250 


500 


35 


250 


517K 


1.72 pb 


300 fb- 1 


2 


300 


500 


35 


300 


494K 


1.65 pb 


300 fb" 1 


3 


350 


500 


35 


350 


470K 


1.57 pb 


300 fb" 1 


4 


400 


500 


35 


400 


442K 


1.48 pb 


300 fb" 1 



Table 5: Parameters of the four CMSSM reference points 1, 2, 3, 4 with increasing values of mo and 
Aq (all mass parameters are given in GeV units). We also quote the numbers of events, the LHC 
cross sections and the assumed sample luminosities. 



Point 




Mu L 


M <i, 




M fl 




Mr 
'if 


Mr 


M h 


1 


1155 


1080 


1084 


388 


242 


206 


314 


422 


116 


2 


1157 


1092 


1095 


388 


282 


206 


355 


453 


115 


3 


1159 


1106 


1109 


388 


323 


206 


398 


487 


115 


4 


1162 


1123 


1126 


388 


363 


206 


442 


524 


115 



Table 6: Relevant sparticle masses (in GeV) for the reference points defined in Table 5. 



largest value of m studied, namely 400 GeV for point 4 above. We conclude that the analysis 
described previously in this paper is quite robust, and the LFV signal has a good likelihood of 
being observable, as long as its branching ratio exceeds about 10%. 
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Figure 10: Observable numbers of LFV /i t t^ pairs for the different points described in Tables 5 and 
6, with increasing values of mo and Aq, for fixed M±/ 2 = 500 GeV and tan/3 = 35. 

6 Conclusions 

Motivated by the neutrino oscillation data, and the ensuing likelihood that low-energy violations 
of charged lepton numbers might be detectable in supersymmetric seesaw models, we have 
explored the observability at the LHC of sparticle decays that violate r lepton number. We have 
worked in the context of the Constrained Minimal Supersymmetric extension of the Standard 
Model (CMSSM) and in a non-minimal SU(5) GUT extension of the theory. Focusing mainly 
on regions of the CMSSM parameters with values of the relic neutralino LSP density that 
fall within the range acceptable to cosmology, we investigate have the SUSY parameter space 
requirements for tau flavour violation to be observable in X2 — ► X + T ± Ai =F decays. We have 
studied the possible signals from hadronic r decays, which we have analyzed at the event 
generator level with the use of PYTHIA. 

Within this framework, we have found the following: 

• The observation of LFV in neutralino decays at the LHC can be possible if T(x 2 — ► Xi T± ^) 
/r(x2 - XirVF) „ o.i. 

• The strong bounds on radiative r decays, as well as the other cosmological and phe- 
nomenological requirements, constrain significantly the allowed parameter space, to the 
extent that the CMSSM and the minimal SU(5) GUT are not promising frameworks for 
observing LFV sparticle decays. 

• Larger ratios can be found in a non-minimal SU(5), where RR slepton mixing may be 
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substantial, enabling the LFV signal to be distinguished clearly from the background. 

• The LFV signal remains observable also at larger values of m than are favoured in the 
usual CMSSM framework. 

We conclude that the search for X2 — > Xi T± ^ decays at the LHC is interesting and comple- 
mentary to the parallel searches for r — > 717 decays, and could be a useful 'canary in the mine' 
for non-minimal GUTs. 
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